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FOREWORD 


The  Earth  Limb  Infrared  Atmospheric  Structure  (ELIAS)  experiment  served  as  a 
test  round  for  the  SPIRIT  program  (SPIRIT  recorded  infrared  measurements  of 
auroral  and  earth-limb  emissions  using  sensors  mounted  on  a  sounding  rocket). 
ELIAS  consisted  of  a  three-channel  scanning  radiometer  aboard  a  high-altitude 
sounding  rocket.  Scientists  launched  the  ELIAS  payload  on  18  March  1983. 

CSE  awarded  Subcontract  No.  85-060  to  Yap  Analytics,  of  Lexington,  MA  obligat¬ 
ing  them  to  develop  a  statistical  analysis  and  characterization  of  the  ELIAS 
experimental  data.  Yap  scientists  completed  their  work  under  this  subcon¬ 
tract  and  submitted  their  final  report  to  CSE  in  December,  1985. 

This  report  presents  and  describes  the  computer  codes  which  Yap  scientists  de¬ 
veloped  to  generate  atmospheric  radiance  profiles  for  various  molecules,  spe¬ 
cies,  and  rotational  bands. 

SUMMARY 

Yap  scientists  successfully  debugged  and  tested  the  infrared  atmosphere  vibra¬ 
tional  temperature  code,  RADT.  They  revised  the  code  to  make  it  flexible  for 
various  applications  and  then  cascaded  it  with  the  Non-Equilibrium  Atmospheric 
code,  NLTE.  The  researchers  then  established  procedures  to  use  the  combined 
codes  to  generate  vibrational  temperature  profiles  and  radiance  profiles  for 
various  molecules,  species,  and  rotational  bands.  Codes  were  than  generated 
to  extract  the  radiance  profiles  from  the  standard  NLTE  program  output  files 
and  transfer  them  from  the  Cyber  system  to  the  Apollo  system.  Finally,  Yap 
created  a  routine  to  plot  the  radiance  profiles'  sum  total  and  SPIRE  program 
measured  data. 
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The  ELIAS  experiment  called  for  launching  a  rocket-borne  complement  of  infra¬ 
red  sensors  into  auroral ly  active  skies  from  Poker  Flat  Research  Ranye  in 
Alaska,  while  ground  observers  simultaneously  recorded  the  auroral  activity  at 
Watson  Lake,  Fort  Nelson,  Peace  River,  and  Fort  Providence,  Canada,  The  prime 

sensor  aboard  the  ELIAS  payload  was  a  telescoped,  liquid-nitrogen-cooled  scan¬ 
ning  radiometer  (NTR-1).  This  radiometer  detected  infrared  radiation  at  9100A 

and  in  the  SWIR  and  MWIR  spectral  regions. 

During  the  ELIAS  mission,  o’.  1  sensors  operated  as  expected  and  the  experiment 

provided  excellent  data.  The  ELIAS  mission  also  verified  the  feasibility  of 
pointing  an  infrared  sensor  manually  with  a  hand-operated  "joystick". 

This  report  describes  the  computer  codes  that  Yap  developed  to  generate  the 
atmospheric  radiance  profiles  for  various  species  and  rotational  bands  of  a 

given  molecule.  This  report  also  describes  the  procedures  which  Yap  estab¬ 
lished  to  streamline  the  various  codes.  The  appendix  contains  listings  of  the 

codes  and  procedures. 


2.  OVERALL  FLOWCHART 


Several  steps  are  involved  in  the  generation  of  the  final  plot  of  the  theoret¬ 
ically  computed  infrared  radiance  versus  altitude  for  the  various  species. 
Figure  1  is  an  example  of  the  final  product.  Included  in  the  plot  are  data 
points  measured  on  the  AFGL  SPIRE  program.  The  radiance  axis  is  in  Watt/Cm2- 
Sr-Band,  and  the  tangent  height  is  in  Km.  Species  data  were  computed  and  il¬ 
lustrated  together  with  the  sum  total  radiances  of  these  species. 

Figure  2  shows  a  general  flowchart  of  the  procedures.  Six  major  procedures 
constitute  the  production  of  the  final  plot.  The  first  four  procedures  are 

performed  on  the  AFGL  Cybe~  system  where  the  major  computer  codes  reside.  A 

transfer  of  the  output  file  is  then  made  from  the  Cyber  to  the  Apollo  system, 

where  the  final  plotting  of  the  band  radiance  profiles  are  then  performed. 

The  first  procedure,  'BOTH',  involves  two  major  sets  of  codes,  as  will  be  de¬ 
scribed  in  detail  in  a  later  section.  The  two  cooes  are  entitled  ' NLTE4 '  and 
RADT ' .  The  main  product  of  the  procedure  'BOTH'  is  a  file  'SPIREF'  contain¬ 
ing  the  vibrational  temperature  versus  altitude  of  a  particular  band  and 
branch  (e.g.  C02,  Q  Dranch).  This  procedure  is  generally  repeated  four  or 

more  times  to  assure  computation  convergence. 

The  vibrational  temperature  profile  is  then  fed  into  a  code  in  the  procedure 
’FIXSPIR1  from  which  the  hot  band  vibrational  temperature  profile  ‘SPIREH1  is 
extracted . 

The  next  procedure,  "JOB#1  ,  is  processed  eleven  times,  once  per  specie  to  ob¬ 
tain  the  infrared  radiance  profile  for  each  specie.  Eleven  procedures  la¬ 

belled  1 J OB 1 '  ,  .......  ' J0B11 '  have  been  set  up  to  perform  the  computations 

for  each  of  the  species. 

Tne  eleven  output  files  from  the  'JOB#'  procedures  are  then  streamed  into  the 
procedure  'PLTFSLE'  which  extracts  just  the  band  radiance  profiles  from  each 
file  and  reformats  them  in  a  concise  file  called  ' X X 0 U T '  for  plotting. 
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At  this  point,  the  reformatted  radiance  profiles  in  file  'XXOUT1  are  trans¬ 
ferred  from  the  Cyber  system  to  the  Apollo  system  through  the  high  rate  N1U 
network,  and  stored  under  the  file  name  ' C02DATA1 ' .  In  the  Apollo,  a  plot 

routine  has  been  formulated  to  accept  the  file  containing  the  radiance  pro¬ 
files,  and  to  plot  each  radiance  profile  with  the  total  radiance,  as  well  as 
the  AFGL  SPIRE  program  measured  data. 


3.  PROCEDURE  ’BOTH 


Procedure  'BOTH'  utilizes  two  major  AFGL  computer  codes  -  NLTE4  (Non- 
Equilibrium  Atmospheric  Code)  and  RADT  (Vibrational  Temperature  Profile  Code). 
Both  codes  are  resident  in  the  AFGL  Cyber  system  and  may  be  accessed  upon  re¬ 
quest. 

Figure  3  shows  the  flowchart  for  Procedure  ’BOTH1.  The  description  and  use  of 

the  NLTE  computer  code  may  be  found  in  the  document  AFGL-TR-83-0168 .  NLTE4  is 
a  modified  version  of  NLTE;  it  provides  a  single  line  profile  as  an  additional 

output  file.  The  input  file  formats  for  Units  1-3  remain  unchanged  from  those 
described  in  the  AFGL  document.  Table  1  is  a  sample  of  Unit  l's  input  file 

format.  In  Card  1C,  the  fifth  and  seventh  parameters  (computation  accuracy 
and  wing  integration)  must  be  set  to  zero  to  obtain  a  complete  output  file  for 
the  single  line  profile. 

The  atmospheric  profile  for  Unit  2  and  the  AFGL  line  profile  for  Unit  3  are 
also  described  in  the  above-mentioned  AFGL  document.  These  may  be  created  by 
following  the  descriptions  in  the  document.  Once  created,  they  must  be  stored 
under  file  names  ' SP I  RE  F '  and  'C0215MI'  respectively,  to  be  accessed  by  the 
NLTE4  program  code. 

The  single  line  profile  at  the  output  of  Unit  7  is  stored  under  the  file  name 
' Z I NPRAD ' .  This  file  is  retrieved  automatically  by  the  ’BOTH’  procedure  to  be 
used  on  Unit  8  as  one  of  the  inputs  to  the  'RADT'  code. 

Program  inputs  to  'RADT'  on  Unit  1  are  illustrated  in  Table  2.  Card  1A  is 
read  in  with  format  ’ A1 *  and  the  two  permitted  line  profiles  are  either  the 

Voight  ('V')  or  the  Doppler  (’O').  Card  IB  is  read  in  with  format  ' A5 '  and 
the  two  permitted  methods  of  computing  the  upwelling  radiance  from  below  are 

either  the  blackbody  ( ’ 8B0DY ' )  or  FASCODE  ( ’ F CODE ' ) .  Card  1C  is  read  in  with 
format  ' A1 '  and  the  permitted  line  branches  are  'P‘,  ' Q ' ,  or  'R'. 


Unit  4  input  to  RADT  is  the  ' SP I  REF '  file  containing  the  atmospheric  profile 
which  is  read  in  with  the  format  described  for  Unit  2  input  to  NLTE4.  This 
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file  is  stored  under  file  name  ’ OLDSP I R 1  upon  completion  of  the  process.  This 

permits  a  comparison  of  the  input  file  * OLDSP I R *  with  the  new  atmospheric  pro¬ 
file  'SPIREF'  at  the  output  on  Unit  7,  to  test  for  computation  convergence. 

The  atmospheric  line  profile  input  on  Unit  5  is  named  'DATAF'.  Table  3  de¬ 
scribes  Unit  5  input  formats.  An  additional  input  record  (la)  is  required  if 
the  choice  for  computing  the  upwelling  radiance  from  below  is  'FCODE1.  In 
this  instance,  a  second  temperature  (T2)  and  a  distance  from  the  line  center 
(N2)  must  be  inserted  as  the  first  input  record. 

An  updated  atmospheric  vibrational  temperature  profile  is  generated  on  Unit  7 
as  output  of  this  procedure.  Its  format  is  identical  to  that  of  input  file 
'SPIREF'  on  Unit  4.  At  the  end  of  the  procedure,  this  file  replaces  the  file 
'SPIREF'.  In  general,  this  entire  procedure  is  cycled  at  least  four  times  to 
assure  convergence  in  the  computation  of  the  vibrational  temperature  profile. 

Test  of  computation  convergence  may  be  verified  by  comparing  the  previous  file 
' OLDSP I R '  with  the  new  file  'SPIREF'.  An  example  of  the  program's  computa¬ 
tional  convergence  is  illustrated  in  Figures  4-5.  In  Figure  4,  the  difference 
between  the  vibrational  temperature  profiles  from  runs  2  and  1  are  depicted  as 

a  precentage  of  the  vibrational  temperature  profiles  from  run  4.  Less  than 
0.03%  difference  was  observed  at  low  altitudes,  while  at  high  altitudes  the 

difference  is  well  below  that.  In  Figure  5  the  difference  obtained  from  runs 
4  and  3  were  again  shown  as  a  percentage  of  that  from  run  4.  Computational 

convergence  over  the  entire  altitude  range  is  evident  in  the  figure. 


•  *• 


The  procedure  is  assembled  as  a  local  file  ' B '  which  is  listed  in  Table  4.  To 
use  procedure  'BOTH',  the  following  instruction  is  entered  on  an  interactive 
terminal:  BEGIN, BOTH, B. 


s»- 
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4.  PROCEDURE  ' FIXSPIR 


The  flowchart  for  procedure  'FIXSPIR'  is  shown  in  Figure  6.  The  vibrational 
temperature  profile  'SPIREF'  obtained  from  cycling  through  procedure  'BOTH' 
(at  least  four  times)  is  placed  on  Unit  10  as  input  to  the  Hot  Band  Vibration¬ 
al  Temperature  computer  code  'FIXSPIR'.  The  code  generates  an  output  file 
'SPIREH'  on  Unit  20.  In  addition  to  a  header  record,  ' SP I  RE  H '  contains  re¬ 
cords  of  the  hot  band  profile,  as  Table  5  illustrates. 

The  format  for  these  records  is  (F5 .0 ,F10.3 ,5E12 .5) .  A  listing  of  the  local 
file  'FIX'  containing  procedure  'FIXSPIR'  is  in  Table  6. 

To  exercise  procedure  'FIXSPIR',  the  following  instruction  is  entered  on  an 
interactive  keyboard:  BEGIN, FIXSPIR,  FIX. 


5.  PROCEDURE  'JOB#' 


Procedure  'JOB#'  (where  #  =  1, . ,11)  utilizes  the  computer  code  ' NLTEE ' 

which  is  an  updated  version  of  Non-Equilibrium  Atmosphere  Code  'NLTE'. 

'NLTEE'  was  revised  to  replace  the  Simpsons'  integration  in  'NLTE'  with 
Gaussian  integration.  Because  of  the  Lorentzian  roll-off  of  the  line  profile, 
the  Gaussian  integration  permits  progressively  coarser  sampling  over  the  wings 
of  the  line  profile,  thus  resulting  in  a  tremendous  time  improvement  in 
running  the  code.  Aside  from  this  improvement,  input  requirements  and  formats 

are  identical  to  those  of  'NLTE'.  The  flowchart  for  the  procedure  is  shown  in 
Figure  7.  One  of  the  following  AFGL  line  profiles  are  used  for  input  on  Unit 

3  for  each  of  the  eleven  successive  runs: 


1. 

' C0215MI ' 

7. 

‘ C02626 ' 

2. 

'C0215MI ' 

8. 

' C02626 ' 

3. 

' C0215MI ' 

9. 

' C02DATA' 

4. 

'  L010000' 

10. 

' C02DATA' 

5. 

' C02626 ' 

11. 

' C02DATA' 

6. 

' C02626 ' 

For  runs  1  through  4,  the  atmospheric  profile  ' SP I REF '  is  used  as  input  on 

Unit  2.  For  the  hot-band  computations 

of  runs  5 

through  11,  the  file  'SPIREH 

is  the  input  on  Unit  2. 

The  band  radiance  profi1 

le  is  outputted 

on 

Unit 

5  with  the  fol lowing 

format:  (F7  .3  ,A3 ,214  ,5(  14  ,E10 .4  ,E9 .3) ) 

fol 

lowing 

1  header  records.  Table  7 

lists  the  output  parameters  on  Unit  5. 

The 

last 

four  parameters  (LBND,  LMAX, 

SUMRD,  TOD)  are  written 

five  times  per 

record . 

The  output  files  from  Unit  5  for 

each 

of  the  eleven  runs  are  stored 

respectively  as: 

1. 

' T562601 * 

7. 

' T562611 ' 

2. 

1 T563601 ' 

8. 

' T562612 ' 

3. 

1 T562801 ' 

9. 

' T56261 1 

4. 

' T562701 ' 

10. 

' T56262 ' 

5. 

' T562603' 

11  . 

' T56302 1 

6. 

' T562602 1 
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Table  8  lists  the  local  file  * B# *  (where  It  =  1,2 . 11)  which  contains  the 

procedure  'JOB#'  (#  =  1,....,11).  To  process  each  of  the  'JOB#'  procedures, 
the  following  command  is  entered  via  an  interactive  keyboard:  BEGIN, JOB# ,B# . 


c* 


6.  PROCEDURE  'PLTFILE 1 

In  the  11-Band  Radiance  Profile  Conversion  Code  1  RE ADFL '  contained  in  the  pro¬ 
cedure  'PLTFILE',  the  eleven  band  outputs  from  the  procedure  'JOB#'  are  con¬ 
verted  into  a  concise  file  suitable  for  plotting.  Each  file,  ('1562601',.. 
..'T563602')  is  sequentially  read  in  on  Unit  10.  Figure  8  is  the  flowchart 
for  the  procedure.  ' READFL '  searches  for  the  band  radiance  for  each  altitude 
from  each  file.  When  all  eleven  files  have  been  processed,  'READFL'  writes 
the  eleven  band-radiance  profiles  out  on  Unit  20.  At  the  termination  of  pro¬ 
cedure  'PLTFLE',  the  radiance  profile  is  stored  under  file  name  ‘ X XOUT '  . 

A  listing  of  the  local  file  'PLTFLE'  containing  the  procedure  'PLTFILE'  is  in 
Table  9.  To  generate  the  concise  radiance  profile  with  procedure  'PLTFILE', 
the  following  command  is  entered  via  keyboard  on  an  interactive  system: 

BEGIN, PLTFILE, PLTFLE. 
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7.  CYBER  TO  APOLLO 

A  plot  routine  has  been  coded  in  Fortran  on  the  Apollo  system  specifically  to 
plot  the  band  radiance  profiles  together  with  their  sum  and  Spire  program 
measured  data.  The  concise  plot  file  has  to  be  transferred  from  the  Cyber  to 
the  Apollo  system  through  the  communication  network  (NIU).  The  plot  file  is 
short  and  takes  a  matter  of  seconds  once  the  link  has  been  established  between 
the  Cyber  and  tne  Apollo.  The  plot  file  on  the  Cyber  ' XXOUT '  is  retrieved  and 
transferred  to  the  Apollo  and  stored  under  the  file  name  ' C02DATA1 *  (see  Fig¬ 
ure  9).  The  file  'C02DATA11  should  be  stored  in  the  following  directory: 

1 /USERS/LSP/BANKS'  since  the  plotting  package  and  necessary  inputs  to  the  plot 
package  are  in  the  same  directory. 


8.  APOLLO  PLOT  ROUTINE 


Figure  10  illustrates  the  plotting  procedure  on  the  Apollo  system.  Unit  8 
contains  the  computed  data  file  transferred  from  the  Cyber.  The  Unit  9  input 
file  contains  the  Spire  program  measurement  data  as  well  as  other  plot  ini¬ 
tial  ization  data. 

The  plotting  program  'co2‘  is  compiled  and  stored  under  file  name  1 BOBEEP ‘ . 
When  1  BOP E EP 1  is  executed,  both  input  files  are  read  and  the  plottiny  coordi¬ 
nates  are  set  up.  The  total  radiance  of  the  eleven  bands  is  computed.  The 
twelve  curves  are  then  plotted  within  the  confines  of  the  setup  coordinates. 

The  plot  routine  on  the  Apollo  is  stored  in  the  directory  of  ‘/USERS/LSP/ 

BANKS'.  To  execute  the  plot  routine,  the  following  instruction  is  entered 
while  working  in  the  above  directory:  ‘BOPEEP1. 

An  output  file  named  1 DATA_PLOT_FILE .01 '  is  generated  after  the  screen  plot¬ 
ting  is  complete.  A  hard  copy  of  this  plot  may  be  obtained  with  the  following 
instructions:  PRF  -NPAG  DATA_PLOT_FILE .01 . 

It  is  advisable  to  delete  the  stored  plot  to  free  up  disk  space  after  the  hard 

copy  is  obtained.  The  following  instruction  is  entered: 

OLF  DATA  PLOT  FILE. 01. 
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F mure  i.  CARBON  DIOXIDE  RADIANCE  PROFILE 
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Table  1  Program  Inputs  on  Unit  1 


Quantity 

Description 

Unit 

Type 

Example 

(Card  1A) 

MOL 

Molecular  Formular 

Character 

'  C02  ’ 

ISO 

Isotope  code 

Integer 

626 

VIBE 

Vib  energy  of  the  transition 

-1 

cm 

Real 

667.379 

VI BL 

Vib  energy  of  lower  state 

—  1 
cm 

Real 

0 

VIBQ 

Vib  energy  for  partition  fn 

-1 

cm 

Real 

667.379 

(Card  IB) 

BR 

Branch  (P,Q,  or  R) 

Character 

•Q' 

NRQL 

Rotational  quantum  line 

Integer 

30 

UST 

Upper  vib  state  (AFGL  not  'n) 

Character 

' 01101' 

LST 

Lower  vib  state  (AFGL  not  'n) 

Character 

'00001' 

NDP 

Lineshape  option 

Integer 

0 

VMIN 

Lowest  wavenumber  desired 

-1 

cm 

Real , int 

600 

VM  AX 

Highest  wavenumber  desired 

cm  * 

Read ,  int 

750 

(Card  1C) 

TANIN 

Initial  tangent  height 

km 

Integer 

65 

TANF 

Final  tangent  height 

km 

Integer 

100 

INTRVL 

Examination  height 

km 

Integer 

5 

HM  AX 

Altitude  of  top  layer 

km 

Integer 

150 

ACC 

Fractional  error 

Rea  1 

0 

NPTS 

Number  of  integration  points 

Rea  1 

WING 

Wing  cutoff 

Real 

0 

Example  of 

:  an  Input  Data  Record  on  Unit 

1 

'C02' ,626, 

667.379,0,667.379/ 

'Q'  ,30,  '01101’ , '00001'/ 
65,100,-1,150,0,  ,0/ 

0 


Table  2.  RADT  Program  Inputs  on  Unit  1. 


Quantity  Description  Unit  Type  Example 


(CARD  1A) 

IANSWER  Line  profile  Character  'V' 

(CARD  IB) 

METHOD  Method  of  radiance  computation  Character  'BBODY' 

(CARD  1C) 

IB  Branch  of  spectral  line  Character  'Q' 


Example  of  an  Input  Data  Record  on  Unit  1  for  RADT 
•  V' 

'BBODY' 


Table  3.  Atmospheric  Line  Profile  on  Unit  5 


Quantity  Description 


(RECORD 

la)  "*****  Required  when  option 

is  FCODE 

T2 

Second  Temperature 

Deg  K 

Real 

N2 

Distance  from  L^ne  Center 

-1 

cm 

Real 

(RECORD 

1) 

AE 

Einstein  Coefficient 

Real 

ROTCON 

Rotational  Constant 

Real 

STWR 

Stat.  Weight  for  Rot.  Line 

Real 

QAL 

Non  zero  if  Q  Branch  allowed 

Real 

STU 

Stat.  Wt .  for  Upper  Vib.  Level 

Real 

STL 

Stat.  Wt .  for  Lower  Vib.  Level 

Real 

(RECORD 

2) 

T 1 

Initial  Temperature 

Deg  K 

Real 

-1 

Real 

DW 

Doppler  Width 

Cm 

S 

Absorption  Coefficient 

Cm/Part. 

Real 

FNU 

Line  Center 

Cm-1 

Real 

HI 

Initial  Altitude 

Km 

Real 

ROT 

Rotational  Quantum  # 

Real 

(RECORD 

3  -  N  ) 

H(N) 

Lower  Altitude 

Km 

Real 

T(  N  ) 

Lower  Temperature 

Deg  K 

Real 

C  (  N  ) 

Lower  CO2  Density 

Part ./CC 

Real 

A  (  N  ) 

Lower  Density 

Part./CC 

Real 

B  (  N  ) 

Lower  0-Atom  Density 

Part./CC 

Rea  1 

CT  (  N  ) 

Lower  0,  Density 

Part./CC 

Rea  1 

Example 


Table  4.  Listing  of  File  'B*  with  Procedure  1 
.PROC,BOTH. 

GET, LG0=NLTE4B 
GET,TAPE1>=NLTE4IN. 

GET,TAPE2=SPIREF. 

GET ,TAPE3=C02 15MI . 

REWIND, TAPE1 , TAPE2 , TAPE3 , OUTPUT, LGO . 

ASSIGN, MS, OUTPUT. 

LGO. 

ASSIGN, TT, OUTPUT. 

REWIND, TXPE7 . 

REPLACE , TAPE7=ZINPRAD. 

RETURN , LGO , TAPEl , TAPE2 , TAPE 3 , TAPE7 , TAPE 5 , TAPE 4 ,TAPE6 
GET , LGO=RADTB . 

GET , TAPEl=RADTIN . 

GET,TAPE4=SPIREF. 

GET , TAPE 5=DATAF . 

GET ,TAPE8=Z INPRAD . 

REWIND, TAPEl, TAPE4 , TAPE5 , TAPE8 , LGO , OUTPUT . 

ASSIGN, MS, OUTPUT. 

LGO. 

ASSIGN, TT, OUTPUT. 

REWIND, TAPE7 . 

REPLACE, TAPE7=SPIREF. 

REPLACE, TAPE4=OLDSPIR. 

RETURN, TAPEl ,TAPE4 , TAPE 5 , TAPE 7 , TAPE 8 , LGO . 


Table  5.  Hot  Band  Profile  'SPIRER 


Quantity  Description 


Example 


ALT 

Altitude 

Km 

Real , int 

TRTMP 

Translational  Temperature 

Deg  K 

Real 

PRESS 

Atmospheric  Pressure 

Atmos . 

Real 

RHO 

#  Density  of  Rad.  Mole. 

Cm-3 

Real 

TVL 

Lower  Level  Vib.  Temp. 

Deg.  K 

Real 

TVU 

Upper  Level  Vib.  Temp. 

Deg  .  K 

Real 

TVQ 

Lowest  Excited  State  Vib. 

Deg .  K 

Real 

Table  6.  Listing  of  Procedure  'FIXSPIR  . 
. P  ROC ,  F I X  S  P I R  . 

GET , F IXS  PB . 

GET , SPI REF . 

REWIND, F IXS PB,S PI REF. 

FIXSPB. 

REWIND, SPIREH . 

REPLACE, SPIREH. 

GET, yY=FIXSPIR. 

RETURN, FIXSPB , SPI REF, SPIREH. 


V  *V  *f  ■*«  •* 
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Table  7.  Radiance  Profile  Output  on  Unit  5. 


Quantity 

Description 

Units 

Type  Sample 

TS 

Real 

BR  (  J ) 

Branch 

Character 

RQL 

Integer 

LBND 

Band 

Integer 

LMAX 

Integer 

SUMRD 

Integrated  Radiance 

W/CM2-Sr 

Real 

TOD 

Table  8.  Listing  of  'B#'  containing  procedure  'JOB#'. 
.PROC. JOB1. 

GET , TAPEl=JOBl IN . 

GET,TAPE2=SPIREF. 

GET/TAPE3=C0215MI . 

GET, LGO=NLTEEB. 

REWIND, LGO,TAPEl,TAPE2 ,TAPE3 , OUTPUT. 

ASSIGN, MS, OUTPUT. 

LGO . 

ASSIGN, TT, OUTPUT. 

REW IND,T APES ,TAPE6 , TAP E7 , OUTPUT. 

REPLACE, TAPE5=T562 601. 

RETURN, TAPE1 ,TAPE2 , TAPE 3 , TAPE 5 , TAPES ,TAPE7 . 

RETURN, LGO. 


TABLE  9.  Listing  of  'PLTFLE'  with  Procedure  'PLTPILE 
.  PROC , PLTFILE . 

GET,XXl=T562601. 

GET,XX2*=T563601. 

GET,XX3=T562801. 

GET,XX4=T562701. 

GET,XX5S=T562603. 

GET, XX6=T562602 . 

GET,XX7=T562611. 

GET , XX8=T562612 . 

GET, XX9=T5626 1 . 

"GET,  XXI  0=T56262 . 

GET,XXll=T563602. 

GET, AA=READB2 . 

REWIND , AA , XXI , XX2 , XX3 , XX4 , XX5  ,  XX6  . 

REWIND, XX7 , XX 8 , XX9 , XX10 , XXll . 

AA. 

RETURN , XXI , XX 2 ,XX3 , XX4 , XX5 , XX6  . 

REWIND, XXOUT. 

RETURN , XX7 , XX8 , XX9 , XX10 , XXll , AA. 

REPLACE, XXOUT. 

RETURN, XXOUT. 


program  co «* 

real  x 626 < 1 8 > > x636<18),  x628<18),  x627<18).  g(18).  total(18) 
real  x 03301 <  18 ) , x02201<18), xlll01<18).  xl 1102(18) 
real  x 10001 ( 18) < xl0002<18>, x63602(18) 
real  sx < 141 ),  sg ( 141 ) 

INTEGER  IC < 18) 
character*10  id(8) 

open  <unit=4, f ile='garbage  ') 
open  ( uni t=8,  f i le« 'co2datal ' ) 
open  <unit=9,  f  i  le=  'data  ' ) 

do  400  n=l » 141 

rea d(9>  2) sg ( n ) , sx(n) 
format(f5.  I,2x,e7.  1) 
continue 


do  200  j  =  l. 18 
total (j)“0.p 
|  continue 

a  1 1=&5.  0 
do  100  i=l,  18 

read <8.  I)x626<i>,  x636(i),  x62B(i),  x627(i),  x03301(i),  x02201(i). 
&  x  1 1 1 0 1  (  i  )  >  x  1 1 1 02 (  i  ) >  x  10001  (i),  xl0002<i>,  x63602(i) 

f  ormatf  1 1  ( 1  x,  elO.  5)) 

total <i >=x626< i )+x636( i )+x628< i )+x627( i >+x03301 ( i )+x02201 ( i )+ 
&  xlll01(l)+xlll02(i)+x 10001 ( i  )  +  x 10002 ( i )+x63602( i ) 

g ( i )=alt 
alt=alt+5.  0 
continue 

DO  10  1  =  1,  18 
IC<I)=1 
10  CONTINUE 

DO  20  1  =  1,  18 

IF  (  X626<  I ) .  GE.  9.  0E-10 )  IC<1)=I 
IF (  X636 ( I ) .  GE.  9.  0E-10 )  IC(2)  =  I 
IF ( X628< I ) .  CE.  9.  0E-10 )  IC(3)=I 
IF( X627 < I ) .  GE.  9.  0E-10 )  IC(4)  =  I 
IF( X03301 ( I ) .  GE.  9.  0E-10)  IC(5)  =  I 
IF  <  X02201 ( I ) .  GE.  9.  0E-10)  IC(6)  =  I 
IF(XUIOKI).  GE.  9.  0E-10)  IC<7)  =  I 
IF(X1U02(  I ).  GE.  9.  OE-IO)  IC(8)  =  I 
IF (X 10001  ( I  ).  CE.  9.  OE-IO)  IC(9)  =  I 
IF ( X 10002 ( I ) .  CE.  9.  OE-IO)  IC<10)  =  I 
IF (  X63602 <  I ) .  GE.  9.  OE-IO)  ZC  C 1 1  >  — I 
IF(TOTAL<I).  CE.  9.  OE-IO)  IC(12)-I 
20  CONTINUE 

call  pltid3<  id,  1,2,  .  7.0) 
cal  1  p  1  ot  < 2.  2,  6.  0.  -3) 


gmin=50.  0 
dg*10.  0 


xmin»-9. 
dx«l.  /3. 
no__dec*5 
zinc-3. 

call  lgaxis(0.  ,  0.  ,  'radiance '.  -B.  15.  .  0.  ,  xmin.  dx,  1.  .  no_dec.  zinc) 
call  lgax is <0.  0*  1 1.  0,  '  '»  1.  15.  0#  0.  0#  xmin.  d x, -1.  #  no_dec#  z inc ) 

call  ax is <0.  0,  0.  0#  'tangent  height '»  14.  1 1.  0*  90.  0,  ymin.  dg<  2.  0) 

call  axis  (IS.  0,  0.  0,  '  ',  -1,  11.  0,  90.  0,  gmin,  0.  0,  2.  0) 

call  symbol  <2.  0,  11.  5,  2.  .  'CARBON  DIOXIDE  RADIANCE  PROFILE'.  0.  31  ) 

call  lgline(sx,  sy >  141.  1.  -1.  15.  xmin.  dx.  ymin.  dy.  2.  1.  0) 
call  lg  1  inc  (  *626,  y<  IC  ( 1  > .  If  1.3.  xmin.  dx.  ymin.  dy.2.  1.0) 

call  lg  1  ine  (  x636,  y>  IC  (2) .  1.  1.2.  xmin. dx, ymin. dy.2.  1.0) 

call  lgline(  x628,  y.  ICO).  1,  1.0.  xmin.dx,  ymin.  dy.2.  1.0) 

call  lgline(  x627.  y,  IC  ( 4) .  1,  1.  5.  xmin,  dx,  ymin,  dy>2>  1,0) 

call  lgline(  x  03301,  y,  IC(5),  1,  1,4,  xmin,  dx.  ymin.  dy.2,  1,0) 

call  lgline<  x  02201,  y,  ICO).  1,  1,8,  xmin,  dx.ymin.  dy.2,  1,0) 

call  lgline(xlll01,  y,  IC  (7) .  1,  1,  13,  xmin.  dx,  ymin,dy,2,  1,0) 

call  lgline(xll  102,  y,  IC  (B ) .  1,  1,  17,  xmin.  dx,  ymin.  dy.2,  1,0) 

call  lgline(  x  10001 ,  y,  IC  ( 9 ) .  1 .  1,6,  xmin,  dx.ymin.  dy»2.  1,0) 

call  lgline(x  10002,  y,  IC(10),  1,  1,7,  xmin.  dx.ymin,  dy.2,  1,0) 

call  lgline(  x 63602,  y,  1C(11),  1,  1.  9,  xmin,  dx,  ymin,  dy.  2.  1,  0) 

call  lgline<total,y,  IC  ( 12 ) ,  1 ,  1,  1, xmin, dx.ymin. dy.2,  1,0) 

call  symbol  <  11.  0.  9.  7,  IB.  ,  'd  '.  0.  0) 

call  symbol  ( 11.  0,  9.  3,  IB.  .  'c  ',  0,  0) 

call  symbol  <  1 1.  0,  8.  9,  18.  ,  'a  0,  0) 

call  symbol  <  11.  0,  8.  5,  18.  .  'f ',0.  0) 

call  symbol  (11.  0.8.  1,  18.  ,  'e',0,  0) 

call  symbol  ( 1 1.  0,  7.  7.  18.  ,  'i  '.  0.  0) 

call  symbol(ll.  0,  7.  3.  18.  ,  'n  0,  0) 

call  symbol  ( 1 1.  O,  6.  9,  18.  ,  'r',0,0) 

call  symbol  ( 11.  0,  5.  3,  18.  ,  'b  '.  0,  0) 

call  symbol(ll.  0,6.  5.  18.  ,  'g',0,0) 

call  symbol  ( 11.  0,  6.  1,  18.  ,  'h  ',  0,  0) 

call  symbol  ( 1 1 .  0.  5.  7,  18.  ,  ' j  0,  0) 

call  symbol  (11.  0.  10.  2,  18.  ,  'p  ',0,  0) 

call  symbol  (11.  5.9.  7.  .  09,  '626  01 101-00001 ',  0,  16) 

call  symbol  ( 1 1.  5,  9.  3,  .  09,  '636  01 101-00001 ',  O,  16) 

call  symbol  (11.  5,8.  9,  .  09,  '628  01 101-00001 ',  0,  16) 

call  symbol  (11.  5,  8.  5,  .  09,  '627  01 101-00001 ',  0,  16) 

call  symbol  (11.  5,8.  1,  .  09,  '626  03301-02201 ',  0,  16) 

call  symbol  (11  5,  7.  7,  .09,  '626  02201-01101 ',  0,  16) 

call  symboldl.  5,  7.  3,  .  09,  '626  1 1 101-10001  ',  0,  16) 

call  symbol  ( 11.  5.  6.  9.  .  09.  '626  1 1102-10002 ',  0,  16) 

call  symbol  ( 11.  3.  6.  5.  .  09,  '626  10001-01 101 ',  0,  16 ) 

call  symbol  ( 11.  5,  6.  1,  .  09,  '626  10002-01 101  ',  0,  16) 

call  symbol  ( 11.  5,  5.  7,  .  09,  '636  02201-01 101 0,  16 > 

call  symboldl.  5.  5.  3,  .  09,  'TOTAL', 0,  5) 

call  symboldl.  5,  10.  2,  .  09,  'SPIRE  DATA '.  0,  10) 

call  endp  It 
end 


PROGRAM  FI XSPIR ( INPUT , OUTPUT ) 
CHARACTER*80  HEADER 
CHARACTER*40  NAME 
NR- 10 
NU-20 

OPEN(UNIT-NR.  FILE- 'SPIREF '  > 

OPEN ( UNI T-NW,  FILE- 'SPIREH' ) 

READ < NR,  80)  HEADER 
80  FORMAT <A80) 

WRITE<NW.  80)  HEADER 
DO  10  1=63, 190 

READ(NR,  '  (A40,  F12.  3)  ' )  NAME,  A 
WRITE(NW,  *  ( A39,  3F12.  3)  ')  NAME,  A,  A.  A 
10  CONTINUE 
CLOSE (NR) 

CLOSE (NW) 

END 


n  r*  n  i 


PROCRAM  RfcADFLC  INPUT.  OUTPUT) 

DIMENSION  AC11.  18).  BO).  CO) 

CHARACTERMOO  D 
INTEGER  KK< 10) 

CHARACTERS  NAME 
DO  110  IJK-1.  11 
IF(IJK.  EO.  1)  NAME*  'XXI' 

IF(IJK.  EQ. 2)  NAM£-'XX2' 

IF(IJK.  EQ.3)  NAME- 'XX3  ' 

IF<IJK.  EQ.  4)  NAME- ' XX4  ' 

IF(IJK.  EQ.  5)  NAME—  ' XX3  ' 

IF(IJK.  EO.  6)  NAME— 'XX6 ' 

IFdJK  EQ.  7)  NAME  -  'XX7  ' 

1F(IJK.  EQ.  8)  NAME-  '  X  X8  ' 

IFdJK.  EQ  9)  NAME-  '  XX9 ' 

IFdJK.  EQ.  10)  NAME—  ' XX 1 0  ' 

IF(IJK  EQ  11)  NAME— '  X  X  1  1 ' 

OPEN ( UN  I T- 1 0,  FILE-NAME) 

DO  50  1-1,11 

READ(  10,  53-  ETRR-30,  END-50)  D 
53  FORMAT ( A 1 00 ) 

50  CONTINUE 

DO  90  J-l,  3 
DO  10  I - 1 , 1000 

READ  (  10,  13,  ERR -20,  END-20)  W.  D.  K.  K.  ( KK  <L>*  B(L),C(L),  L—  1 «  3 ) 
13  FORMAT  (F7.  2,  A2-  I  3,  1 4.  3  (1 4,  £  1 0.  4,  E9  3)  ) 

IF  < W  EQ  0  0)  COTO  20 

DO  25  L«=l,  3 

Ad  JK.  J*5-3+L  )=8  (L  ) 

25  CONTINUE 
10  CONTINUE 
20  CONTINUE 

DO  28  L«=l,  6 

READ(  10,  35.  ERR-28,  END-28)  D 
28  CONTINUE 
90  CONTINUE 

DO  40  1=1 , 1000 

READ  (  10,  4  5,  ERR  =  60,  END=60>  W,  D,  K,  K,  ( KK  ( L  > ,  B  (  L  ) .  C  ( L  ) ,  L- 1 . 3 ) 
4  5  FORMAT  <F7.  2,  A2,  I  5,  1 4,  3  (1 4 ,  E  1  0.  4,  E9  3)  > 

IF(W  EQ  0  0)  COTO  60 
DO  75  L- 1 , 3 

A ( I JK,  1 5  +  L ) -B ( L  ) 

75  CONTINUE 
40  CONTINUE 
60  CONTINUE 
CLOSE ( 10) 

110  CONTINUE 

OPEN ( UN  I T-20,  FILE- 'XXOUT '  ) 

DO  70  1-1-18 

WR  I  TE  <  •  -  '  (  1  1E1  1  3)  '  )  <  Ad  JK.  I),  IJK-1.  11) 

WR I TE <20.  '<  2  IE  1  1  5)')  ( A d JK .  I ) ,  I JK- 1 .  1 1  ) 

70  CONTINUE 
CLOSE  <  20 ) 

CND 


